Gap Anisotropy in Iron-Based Superconductors: A Point-Contact Andreev Reflection 

Study of BaFe 2 _ x .Nl e As 2 Single Crystals 
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We report a systematic investigation on c-axis point-contact Andreev reflection (PCAR) in 
BaFe2- 2! Ni a; As2 superconducting single crystals from underdoped to overdoped regions (0.075 
< x < 0.15). At optimal doping (x = 0.1) the PCAR spectrum feature the structures of two 
superconducting gap and electron-boson coupling mode. In the s± scenario, quantitative analysis 
using a generalized Blonder- Tinkham-Klapwijk (BTK) formalism with two gaps: one isotropic and 
another angle dependent, suggest a nodeless state in strong-coupling limit with gap minima on the 
Fermi surfaces. Upon crossing above the optimal doping (x > 0.1), the PCAR spectrum show an 
in-gap sharp narrow peak at low bias, in contrast to the case of underdoped samples (x < 0.1), 
signaling the onset of deepened gap minima or nodes in the superconducting gap. This result pro- 
vides evidence of the modulation of the gap amplitude with doping concentration, consistent with 
the calculations for the orbital dependent pair interaction mediated by the antiferromagnetic spin 
fluctuations. 

PACS numbers: 74.20.Rp, 74.25.Ha, 74.70.Dd 



It is generally accepted that superconductivity in iron 
pnictides results from a superexchange repulsion me- 
diated by magnetic excitations, which couple electron 
and hole pockets of the Fermi surface [H-IH- Such pair- 
ing interactions favor either isotropic s-wave order pa- 
rameters with opposite signs on different sheets of the 
Fermi surface (FS)(s± model) or anisotropic s-wave or 
even d-wave order parameters with nodes 0, @j. Con- 
sensuses have been reached on several systems, e. g. 
LaFePO @, KFe 2 As 2 0, BaFe 2 (Asi_ a; P :c )2 8], and so 
on, that nodes exist on the gap structure. However, ex- 
perimental confirmations of such a nodal-gap state re- 
mains highly controversial in other systems }9H20l|. For 
example, measurements of the electronic specific heat of 
Ba(Fei_ :E Co a ;)2As2 have shown a field dep endence con- 
sistent with both a fully gapped FS 12| and a nodal 
quasiparticals at the Fermi level 2 12 1- Such scattered 
experimental results and interpretations may come from 
the different qualities and doping level of the samples 
studied. 

Point-contact Andreev reflection (PCAR) spec- 
troscopy has been adopted for probing the density of 
state (DOS) of superconductors with the high energy res- 
olution. In addition, the capability of this technique to 
study the anisotropy and the temperature dependence 
of the superconducting gap make it a unique tool in 
providing invaluable information for various mechanisms 
of unconventional superconductivity (for a review, see 
Refs. 21, 22). Several theoretical calculations have been 
reported on the PCAR conductance characteristics of 
a junction involving the s± symmetry in iron pnictide 
superconductors [23|, [24 \. However, due to the long- 




standing issue of surface or/and interface degradation, 
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FIG. 1: (Color online) Phase diagram of BaFe2- :E Ni a ;As2 as 
a function of Ni concentration x. The orthorhombic phase 
below T s and the antiferromagnetic (AF) phase below Tn are 
also shown here. The arrows indicate the doping levels of the 
samples under investigation. Inset: Temperature dependence 
of the in-plane resistivity p for samples with the Ni nomi- 
nal doping level x as labeled. Data are vertically shifted for 
clarity. 



experimental results by PCAR technique reveal a wide 
variation in the measured Andreev conductance spectra 
and consequently, the gap values, especially for the case 
of c-axis junctions [25l |26|. In this Letter, we fabricate 
highly transparent c-axis direct contacts to perform the 
PCAR spectroscopy study on a series of electron-doped 
BaFe2_ x Ni 2 ,As 2 single crystals over a wide doping range. 
The conductance spectra show a systematic and consis- 
tent behavior with the variation of the doping level, in- 
dicative of a doping dependence of the order parameter 
for Ni-122 superconductor. In the s± scenario, by using a 
generalized two-gap Blonder-Tinkham-Klapwijk (BTK) 
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FIG. 2: Normalized conductance curves at T = 2 K for c-axis 
contacts (a) OP20a and (b) OP20b. The red solid lines are 
their two-gap fits with the relevant fitting parameters. The 
blue arrows indicate the additional conductance peak at the 
edge of gap and the dashed lines mark the "kink" structures 
in the main-gap G(V) curves. Insets: the raw G(V) curves 
in temperatures from 2 K to 21 K in a step of 2 K. Data 
are vertically shifted for clarity. The black arrows mark the 
corresponding edge-gap conductance peaks in the raw G(V) 
curves for OP20a and OP20b, respectively. 



model, we estimate the gap amplitude on the hole and 
electron FS sheets. 

High-quality single crystals of BaFe2- :c Ni a ;As2 were 
gro wn from a self-flux method, as described elsewhere 
27j . The crystals were characterized using x-ray diffrac- 
tion and energy dispersion (EDX). The doping level in 
the crystals was determined by inductive coupled plasma 
emission spectrometer (ICP), which gave a Ni concentra- 
tion roughly 0.8 times the nominal content x. We choose 
five compositions: underdoped, cc=0.075 (UD13), 0.085 
(UD17); overdoped, with a;=0.12 (OD18), 0.15 (OD14), 
and optimally doped with x = 0.1 (OP20). The typ- 
ical level of impurity phases has been checked by spe- 
cific heat measurement on the optimally doped crystal 
x = 0.1, in which a residual component 70 at T — > 
revealed an impurity phases of ~ 4% 28]. The tempera- 
ture dependence of resistivity for these five compositions 
under investigation is displayed in inset of Fig. 1, by 
which the bulk transition temperature T c is determined 
(95% of the normal state resistivity) for each composi- 
tion. Consequently, the T c value for each composition is 
shown on the phase diagram in the main panel of Fig. 1. 

Point contacts to the flat and shiny surfaces cleaved 
along the c-axis of BaFe2- a; Ni :I ;As2 crystals were made us- 
ing thick silver paste (4929N DuPont) bonding with gold 
wires (of 16 /jm diameter). The typical size of these pla- 
nar contact is about 0.08-0.15 mm under a microscope. 
Due to the nanocrystalline nature of the silver paint, the 



contact made in this way, is actually formed by many 
nanocontacts analog to tip point-contact technique [22j . 
For the backside electrical wiring, we applied ultrapure 
indium or silver paste to cover the whole area of the bot- 
tom surfaces of the crystals. On each piece of the crystal, 
5-6 planar contacts were made from point to point to en- 
sure the reproducibility and consistency of the junction 
conductance spectra and their spectroscopic nature. 

Fig. 2(a) and (b) show the raw (inset) and normal- 
ized conductance curves G(V) — dI(V)/dV of two c-axis 
Ag/BaFei.gNi .iAs 2 point contacts (OP20a and OP20b), 
respectively. The contacts made in this way remain sta- 
ble in thermal cycling, and the contact resistance at high 
bias Rn varies very little (< 6%) over the whole T range 
up to T c . The Andreev signal as the conductance en- 
hancement decreases on increasing T and vanishes at 
T > T c , leaving a slightly asymmetrical ^-shaped normal 
state. Shown in the main panels of Fig. 2(a) and (b), 
the magnitude of the Andreev reflection reaches as high 
as 30%-40%, implying a relatively transparent bound- 
ary between Ag nanoparticle and BaFe2- ;E Nl ! ;As2 super- 
conductors. The stabilities in i?/v and the high level 
of Andreev signal indicate that the conduction channels 
through the contact is in ballistic regime, and therefore, 
energy-resolved spectroscopy is possible. A feature shows 
up in these conductance curves: an additional peak at 
~ 20 mV, and the peak gradually disappears with T 
approaching T c . It seems that this peak is much pro- 
nounced when the Andreev signal is relatively low, which 
is close to the case of tunneling side (see below) . Very re- 
cently, this conductance peak at the edge of the gap has 
been observed in Co-122 crystals [29j], and is attributed 
to the signature of an electron-boson coupling associated 
with the superconducting gap. The observation of the 
electron-boson coupling mode in the conductance spec- 
tra implies the high quality of the point-contacts and thus 
their spectroscopic nature. 

The two-gap superconductivity manifests itself as a 
"kink" in the in- gap conductance, marked at the dashed 
lines in Fig. 2 (a) and (b). To explicitly describe 
the variety of spectral behavior observed and quanti- 
tatively resolve the gap amplitude, we invoke a gener- 
alized BTK formula [30j with three parameters: a di- 
mensionless parameter Z which represents the interface 
transparency; an imaginary quasiparticle energy modifi- 
cation r 31[ which reflects the spectral broadening, and 
the superconducting gap A. In BTK model, the nor- 
mal and Andreev reflection probabilities, respectively, 
are related to the DOS of the superconductor N s — 
NnRe(—. S^M ) with N the normal-state DOS 

and 9 the crystalline angle parallel to the current injec- 
tion. To choose a gap function to calculate these two- 
gap conductance spectra, we assume, based on the s± 
scenario, an isotropic gap and an anisotropic gap 
of the general form A e [l — r + rcos(20)], with the gap 
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FIG. 3: (Color online) Temperature dependence of the nor- 
malized conductance spectra and the relevant two-gap BTK 
fits (solid red lines) for (a) OP20a, and (b) OP20b, respec- 
tively. Data and their fits are vertically shifted for clarity 
except the bottom ones. The step of T increase is 2 K. The 
obtained gap magnitude for (b) OP20a, and (d) OP20b, re- 
spectively, as a function of T. The solid lines are the fits to 
an empirical gap function, see text. 



(isotropic s± state) 



32, 33 1 . Therefore, 



anisotropy ratio r varying from r = 
to r = 1 (completely ci-wave) [2C 
by the standard two component conductance (current) 
model the conductance spectra is the contributions of 
the hole- like (Gh) and electron-like (G e ) Fermi pockets: 
G = wGh + (1 — w)G e , where w is the spectral weight. 
For simplicity, we assume a balanced contribution of hole 
and electron Fermi surfaces to the spectral conductance 
by taking w—0.5. 

Examples of normalized G(V) curves and their fits at 
T = 2 K are shown in the main panels of Fig. 2(a) 
and (b) for junctions OP20a and OP20b, respectively. 
The two-gap BTK model (red lines) fit very well the 
main features of the experimental G(V) curves except 
the electron-boson coupling mode around 20 meV, yield- 
ing a set of fitting parameters associated with gap mag- 
nitude and anisotropy ratio: A/, = 10.5 meV, A e =4.2, 
and r=0.3 for OP20a and A h = 9.5 meV, A e = 4.5 
meV, and r=0.3 for OP20b, respectively. It is noted that 
the same gap magnitude is also extracted from a recent 
PCAR experiment on a c-axis Ag/BaFei.sCoo.2As2 with 
the comparable T c = 24 K [2j|. Here we emphasize that 
a two-gap formula with two s-wave gap (r = 0) can also 
fits our experimental data rather well. However, the pa- 
rameters Y/ A — 0.6 — 0.7 are applied to fulfill the fit at 
this low T, which brings a large uncertainty in the gap 
magnitude. 

With these fitting parameters, we check the validity of 
these fits by extending the fit to the overall temperature 
spectral. As shown in Fig. 3(a) and (c), the two-gap s± 
model still fits reasonably well the T-dependence of these 
G(V) curves with fitted gap magnitude. In this overall- 



and Y h = 2.9 and T e = 0.8 meV (OP20b) are almost 
constant or slightly increase with T. From the fits of 
various curves we obtain the gaps A^ and A e as a func- 
tion of T, which is plotted in Fig. 3(b) and (d) for these 
two junctions, respectively. For comparison, the obtained 
gaps can be approximated by an empirical gap formula: 
A(T) = A tanh(a^/T c /T - 1) with a = 1.95 for A h and 
1.86 for A e (cf. a = 1.74 for weak-coupling BCS gap). 

We analyze the physical meanings of the obtained 
gap values and gap function. It is shown from angle- 
resolved photoemission spectroscopy experiment on a Co- 
122 crystal that the large gap Ah is located on the hole FS 
sheet, instead, the small gap A e is presented on one of the 
electron FS sheets [l(|- The gap values 2A/,/fcgr c s» 11.6 
and 2A e / ksT c w 5.0, both above the BCS weak-coupling 
ratio. Besides, the a value from the A(T) function also 
points to a strong-coupling character for both (hole 
FS) and A e (outer electron FS). These results are consis- 
tent with a three-band s± Eliashberg model [3] , in which 
spin fluctuations mainly provide the interband coupling, 
and thus so in the electron-boson coupling matrix. On 
the other aspect, the existence of strong electron-boson 
coupling in this compound is manifested by the obser- 
vation of the spectral peak E p at about 20 mV. In our 
low-transparency (large Z = 0.3 for OP20a) point con- 
tact, a characteristic energy of Q;, = E p — A max = 13 
meV and 11 meV (OP20b junction with small Z = 0.2). 
This energy scale is compatible with the spin-resonance 
energy observed by neutron scattering on the same crys- 
tals (351 ] . 

The obtained anisotropy ratio r = 0.3, resolved in our 
c-axis PCAR spectroscopy of BaFe1.9Nio.1As2, indicates 
a full gap state with gap minima along c axis. This node- 
less state of optimally doped Ni-122 is in similarity with 
that of optimally-doped Co-122, in which a gap minima 
is already present at maximal T c by the c-axis thermal 
conductivity measurements [lH HH . 

It is natural, however, to inspect the manner of the 
superconducting gaps in the crystals with doping away 
from the optimum. We have measured the point-contact 
G(V) curves in whole T range up to T c for junctions with 
x = 0.075 (UD13), 0.085 (UD17), 0.12 (OD18), and 0.15 
(OD14). The typical G(V) curves at T = 2 K (< 0.2 T c ) 
and T > T c are shown in Fig. 2(a)-(d) for these four 
samples respectively. As shown, these G(V) curves ex- 
hibit a consistent behavior: 1) An underlying feature of 
a dominant single gap is unambiguously identified with 
a similar conductance enhancement of 25%-35% for each 
junction; 2) A parabolic normal-state G(V) curve with 
a slight asymmetry at T > T c for each x, opposite to 
those of hole-doped K-122 
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T spectral fit, r = 0.3 and Zhi e \ - 
stant with T while Th = 3.3 and T P 



0.3 — 0.2 are con- 
= 1.3 meV (OP20a) 



implies the similar 
origin of the underlying normal-state background. Nev- 
ertheless, a striking feature in these normalized (and the 
raw) G(V) curves is that: at T = 2 K, a conductance 
plateau and/or a double peak around zero bias for junc- 
tions UD13 and UD17 gradually evolutes into an in-gap 
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FIG. 4: Normalized conductance curves at T — 2K for (a) 
UD13, (b) UD17, (c) OD18, and (d) OD14. Insets: The 
corresponding raw data of conductance spectra at 2 K and 
T>T C . 



from the analytical fitting, the small gap on the electron- 
like FS sheets shows a crossover from a nodeless in the 
underdoped side to a nodal feature in the overdoped re- 
gion. This result provides evidence of the modulation of 
the gap amplitude on the FS with doping concentration, 
consistent with the calculation for the orbital dependent 
pair interaction mediated by the antiferromagnetic spin 
fluctuations. 
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sharp peak in G(V) for junctions OD18 and OD14. Con- 
sidering the overall spectral consistency in these junc- 
tions, the systematic evolution of the Andreev conduc- 
tance spectra with doping concentration is nontrivial. 
Qualitatively, for highly transparent junctions at finite 
T , the appearance of an in-gap plateau in Andreev con- 
ductance spectrum is a signature of a fully gapped state. 
In contrast, an in-gap conductance peak is a character- 
istic of an anisotropic gap state due to the presence of a 
finite DOS at low energy, like a <i-wave gap in cuprates 
21j . We note that it is not easy to describe the spectral 
behavior using simple formulism, because we are dealing 
with multi-band or/and even multi-gap system. There- 
fore, our observation that the systematic evolution from 
the in-gap conductance plateau for the underdoped sam- 
ples to the in-gap peak in G(V) curves for the overdoped 
samples indicates the existence of doping induced evo- 
lution of superconducting gaps with an isotropic feature 
in the underdoped region to an anisotropic, even, nodal 
gap in the overdoped side. This is highly consistent with 
the result of the T-dependent penetration depth A in a 
series of Ni-122 superconductors, in which AA oc T" with 
the exponent n > 2 for underdoped samples and AA be- 
comes more linear-T dependent for overdoped samples, 
indicating the development of nodal gaps in the over- 
doped region 

In conclusion, measurements on point-contact junc- 
tions made on single crystals of BaFea;Nii_ x As2 illustrate 
an interesting evolution of the gap structure. The An- 
dreev conductance spectra clearly show a full-gap state 
for underdoped crystals and a highly anisotropic, per- 
haps nodal-like gap state for overdoped crystals. Quan- 
titative analysis of the spectral data of optimally-doped 
contacts using a generalized BTK formalism resolves two 
superconducting gaps in strong coupling limit. Resulted 
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